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Graph Induced Complex on Point Data 

Tamal K. Dey* Fengtao Fan* Yusu Wang* 



Abstract 

The efficiency of extracting topological information from point data depends largely on the complex 
that is built on top of the data points. From a computational viewpoint, the most favored complexes for 
this purpose have so far been Vietoris-Rips and witness complexes. While the Vietoris-Rips complex 
is simple to compute and is a good vehicle for extracting topology of sampled spaces, its size is huge- 
particularly in high dimensions. The witness complex on the other hand enjoys a smaller size because of 
a subsampling, but fails to capture the topology in high dimensions unless imposed with extra structures. 
We investigate a complex called the graph induced complex that, to some extent, enjoys the advantages 
of both. It works on a subsample but still retains the power of capturing the topology as the Vietoris-Rips 
complex. It only needs a graph connecting the original sample points from which it builds a complex on 
the subsample thus taming the size considerably. We show that, using the graph induced complex one 
can (i) infer the one dimensional homology of a manifold from a very lean subsample, (ii) reconstruct 
a surface in three dimension from a sparse subsample without computing Delaunay triangulations, (iii) 
infer the persistent homology groups of compact sets from a sufficiently dense sample. We provide 
experimental evidences in support of our theory. 
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1 Introduction 

Acquiring knowledge about a sampled space from a point data has become a key problem in many areas of 
science and engineering. The sampled space could be a hidden manifold sitting in some high dimensions, 
or could be a compact subset of some Euclidean space. Topological information such as the rank of the 
homology groups, or their persistent behavior can divulge important features of the hidden space. Therefore, 
a considerable effort has ensued to extract topological information from point data in recent years ||7] |9] [14] 
[T8l . With the advent of advanced technologies, the data is often generated in abundance. Mixed with the 
burden of high dimensionality, large data sets pose a challenge to the processing resource. As a result, some 
recent investigations have focused on how to use a lighter data structure or sparsify the input which aids a 
faster computation, but still guarantees that the output inference is correct. 

Point data by themselves do not have interesting topology. So, a foremost step in topology inference is 
to impose a structure such as a simplicial complex onto it. The Delaunay, Cech, Vietoris-Rips, and witness 
complexes are some of the most commonly proposed complexes for this purpose. Among these, Vietoris- 
Rips (Rips in short) and witness complexes (9J have been favored because they can be constructed with 
simple computations. Rips complexes are easy to construct as they can be built from a graph by recognizing 
the cliques in it. However, the presence of simplices corresponding to all cliques makes its size quite large. 
Even in three dimensions with a few thousand points, the size of the Rips complex can be an obstacle, if not 
a stopper, for further processing. Witness complexes, on the other hand, have too few simplices to capture 
the topology of the sampled space in dimensions three or more EJ. To tackle this issue, Boissonnat et 
al. H suggested modifications to the original definition of witness complex [23]. This enlarges the witness 
complexes but makes it more complicated and costly to compute. 




Figure 1: A graph induced complex shown with bold vertices, edges, and a shaded triangle on left. Input 
graph within the shaded triangle is shown on right. 



We investigate a new complex, a version of which was originally introduced in |[T6l for the application 
of sensor network routing. We set up a more general definition and call it the graph induced complex. We 
provide new theoretical understanding of the graph induced complex in terms of topology inference. In 
particular, we show that, when equipped with appropriate metric, this complex can decipher the topology 
from data. It retains the simplicity of the Rips complex as well as the sparsity of the witness complex. Its 
construction resembles the sparsified Rips complex proposed in ll22l and also the combinatorial Delaunay 
triangulation proposed in ||6], but it does not build a Rips complex on the subsample and thus is sparser 
than the Rips complex with the same set of vertices. This fact makes a real difference in practice as our 
preliminary experiments show. The idea of graph induced complex also bears similarity to the geodesic 
Delaunay triangulation which was proposed to recover the topology of a bounded planar region (with holes) 
from point samples ifTTl . Our work extends it, as well as investigates its theoretical properties, to more 



general setting beyond the planar case. 

Given a graph G on a point data P equipped with a metric, one can build a graph induced complex 
on a subsample Q C P by throwing in a simplex with a vertex set V C Q if a set of points in P, each 
being closest to exactly one vertex in V, forms a clique. Figure Q] shows a graph induced complex for a 
point data in the plane. Subsampled points are the darker vertices. Input points are grouped according to 
the proximity to the subsampled vertices (indicated with a Voronoi partition). The shaded triangle enlarged 
on right is in the graph induced complex since there is a 3-clique in the input graph whose 3 vertices have 3 
different closest point in the subsample. Observe that, in this example, the graph induced complex has the 
same homology as the sampled space. 

Figure [2] shows experimental results on two data sets, 40,000 sample points from a Klein bottle in M 4 
and 15,000 sample points from the primary circle of natural image data considered in M 25 U. The graphs 
connecting any two points within a = 0.05 unit distance for Klein bottle and a = 0.6 unit distance for the 
primary circle were taken as input for the graph induced complexes. The 2-skeleton of the Rips complexes 
for these a parameters have 608, 200 and 1, 329, 672, 867 simplices respectively. These sizes are too large 
to cany out fast computations. 
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Figure 2: Comparison results for Klein bottle in R 4 (top row) and primaiy circle in M 25 (bottom row). The 
estimated /3i for three complexes are shown on the left, and their sizes are shown on log scale on right. 



For comparisons, we constructed the graph induced complex, sparsified Rips complex, and the witness 
complex on the same subsample determined by a parameter 5. The parameter 5 is also used in the graph 
induced complex (see definitions later) and the witness complex. The edges in the Rips complex built on the 
same subsample were of lengths at most a + 25 (justified by Proposition I2.8I ). We varied 5 and observed the 



rank of the one dimensional homology group (/3i). As evident from the plots, the graph induced complex 
captured fi\ correctly for a significantly wider range of 5 (left plots) while its size remained comparable to 
that of the witness complex (right plots). In some cases, the graph induced complex could capture the correct 
j3i with remarkably small number of simplices. For example, it had $\ = 2 for Klein bottle when there were 
278 simplices for 5 = 0.7 and 154 simplices for 5 = 1.0. In both cases Rips and witness complexes had 
wrong /3i while the Rips complex had a much larger size (log e scale plot) and the witness complex had 
comparable size. This illustrates why the graph induced complex can be a better choice than the Rips and 
witness complexes. 

We establish three different results. First, we show that the one-dimensional homology group of sur- 
faces in three dimensions can be determined by graph induced complexes. Even the surface itself can be 
reconstructed with some post-processing from a sparse subs ample of a sample that could be excessively 
dense. Second, we show that, for higher-dimensional manifolds, one-dimensional homology can still be 
determined from graph induced complexes with a simple modification of the metric. Finally, we extend our 
results to other homology groups where we show that the persistent homology groups of a pair of graph 
induced complexes can determine the homology groups of compact spaces. Experimental results support 
our theory. 

2 Graph induced complex and preliminaries 

First we define the graph induced complex in a more abstract setting which does not require any metric. 

Definition 2.1 Let G(V) be a graph with the vertex set V and let v : V — )■ V be a vertex map where v(V) = 
V C V. The graph induced complex Q(V, V , v) is defined as the simplicial complex where a k-simplex 
a = {v[, v' 2 , . . . , v' k+1 } is in GiV, V, v) if and only if there exists a (k + l)-clique {v\,V2, • • • , ffc+i} Q V 
so that v{vi) = v'ifor each i € {1,2, . . . ,k + 1}. To see that it is indeed a simplicial complex, observe that 
a subset of a clique is also a clique. 

Now we specialize the graph induced complex to the case where the vertices of the input graph comes 
from a metric space. 

Definition 2.2 A metric space (X, d) is a tuple where X is a set and d : X x X — > R + is a distance function 
satisfying d(x, y) > 0, d(x, y) = iff x = y, d(x, y) = d(y, x), and d{x, y) < d(x, z) + d(z, y). 

Definition 2.3 Let (P, d) be a metric space where P is a finite point set and let Q C P be a subset. Let 
Vd : P — )> Q denote the nearest point map where Vd{p) is a point in argmin g q d(p, Q). Given a graph 
G(P) with P as its vertex set, we define its graph induced complex as Q{P, Q, d) := Q(P, Q, I'd)- 

Among the many possible choices for d, we will focus on two cases where d = ds, the Euclidean 
distance, and d = do, the graph distance induced by the graph G(P) assuming its edges have non-negative 
weights. For any two vertices p\,p2 G P, the distance dc(pi,P2) is the length of the shortest path between 
pi and p2 in G(P). We will describe the choices of the distance functions as and when necessary. 

In our case, the point set P will be a discrete subset of a compact smooth manifold M C W 1 without 
boundary, or simply of a compact set X C W l . The graph G(P) will be the graph G a (P) = (P, E a ) where 
(PI1P2) G E a if and only if ||pi — P2II < en. The graph induced complex induced by G a (P) on a subset 
Q C P under a distance function d will be the focus of our study. To emphasize the dependence on the 
parameter a, we denote it as Q a {P, Q, d). One may draw a parallel between the graph induced complexes 
and the well-known witness complexes [231 where P acts as a witness set W and Q acts as a landmark set 
L C W . However, the analogy does not extend any further since the construction of the witness complex 
and its variants H differs from that of the graph induced complex. For example, the original witness 



complex defined in E3l embraces a /c-simplex with vertex set in L only if its vertices belong to the k- 
nearest neighbors of a point in W. In contrast the graph induced complex embraces a fc-simplex only if 
its vertices have nearest neighbors in W that form a fc-clique in a graph built on the vertices belonging to 
W. Similar to the witness complexes, the graph induced complex builds upon a subsampling. But, unlike 
witness complexes, it enjoys some topological properties without any extra modifications such as weighting 
0. 

2.1 Sampling, homology, and sandwiching 

As indicated before, the input point set P is a sample of a smooth manifold M or of a compact set X 
embedded in an Euclidean space. We will also subsample P according to a distance function d. Therefore, 
we define sampling in a more general context. 

Definition 2.4 A finite set P C X is an e-sample of a metric space (X, d), if for each point x £ X there 
is a point in p € P so that d(x, p) < e. Additionally, P is called 5-sparse if d{p\ , P2) > 5 for any pair of 
points in P. 

The point set P does not have interesting topology by itself. We build simplicial complexes using P as 
the vertex set to infer the topology of the sampled space X. Specifically, our goal is to infer the homology 
groups of a manifold or a compact set from which P is sampled by computing the homology groups of a 
simplicial complex built with P as vertices. Let H r (-) denote the r-dimensional homology group. It refers 
to the singular homology when the argument is a manifold or a compact set, and to the simplicial homology 
when it is a simplicial complex. Also, all homology groups are assumed to be denned over the finite field 
Z 2 . 

Our main tool for topological inference rests on the relations of the graph induced complexes to the Rips 
complexes that are known to capture information about the homology groups of spaces 



Definition 2.5 Given a point set P C IR n and a parameter a, the Rips complex lZ a (P) = 1Z a (P, c^) is a 
simplicial complex where a simplex a € lZ a {P) if and only if all vertices of a, drawn from P, are within a 
Euclidean distance of each other. 

Notice that we define Rips complexes with Euclidean distances instead of general metrics which will be 
assumed throughout this paper. It is known that such Rips complexes capture the topology of a manifold M 
if the parameters are chosen right ||2] |20]. We utilize this fact to infer V\\{M) by exploiting a sandwiching 
property of graph induced complexes by Rips complexes. To prove this fact, we recall the concept of 
contiguous maps from algebraic topology. Our main interest in this concept is the fact that two contiguous 
maps between two simplicial complexes induce the same homomoiphism at the homology level. 

Definition 2.6 ( H21I ) Let K,\ and IC2 be two simplicial complexes connected by two simplicial maps a : 
K,\ — > K2 and b : K,\ — > YZ2. We say a and b are contiguous, if and only if for any simplex a € K,\, the 
simplices a(a) and b(a) are faces of a common simplex in /C2. 

Fact 2.7 ( H21I ) If a : K,\ — > /C2 and b : K\ — >• IC2 are contiguous, then the induced homomorphisms 
a* : H r (/Ci) — > H r {K-2) and b* : H r (/Ci) — > H r (/C2) are equal. 

In our case two simplicial complexes will be K,\ = lZ a (P) and /C2 = TZ^(P) for some (3 > a. The map 
a is an inclusion lZ a (P) ^-> 1Z"{P). For the map b, we consider a simplicial map h : lZ a {P) — > G a (P, Q, d) 
which composed with an inclusion Q a (P, Q, d) ■— > TZ^{P) provides b. We elaborate on this construction. 

The vertex sets of 1Z a (P) and G a (P, Q, d) are P and Q respectively with a vertex map v : P — > Q 
where p G P maps to one of its closest point v(p) € Q with respect to the distance function d. Observe 



that G a {P) is the 1-skeleton of lZ a (P). Therefore, the edges of a (k + l)-clique in G a (P) constitute 
the 1-skeleton of a fc-simplex in lZ a (P) and vice versa. The vertex map v extends to a simplicial map 
h : 7Z a (P) — > Q a (P,Q,d) where a fc-simplex {pi,P2, • • • ,Pk+i} in TZ a {P) is mapped to a simplex (of 
dimension at most k) with the vertex set UilKPi)}- To see that h is well defined, observe that any subset 
of the (k + 1) -clique {pi,P2, ■ ■ ■ ,Pk+i} is also a clique in G a (P) and hence UjIKk)} is a simplex in 
Q a (P, Q, d). The following result is used later. 

Proposition 2.8 Let (P, d) be a metric space where P C M. n is a finite set and for every pair pi,P2 £ P, 
d(pi , P2) is at least the Euclidean distance \\p\ — P2 1|. Let Q be a 5-sample of (P, d). We have the sequence 

K a {P) A G a (P,Q,d) A K a+2S (P) 

where j is an inclusion and j oh is contiguous to the inclusion i : lZ a (P) ^-> lZ a+ (P). Hence, j*oh* = i*. 

Proof: The map h is well-defined as we detailed before. We observe that Q a (P, Q, d) C lZ a+2S (P) because 
any edge (91,92) of a simplex a G Q a (P,Q,d) satisfies d(qi, 92) < a + 25. Since ^(91,92) > ||9i — 92II by 
assumption, the edge (91, 92) and hence the simplex a wee in lZ a+2S (P). It follows that the inclusion map j 
is well-defined. 

To prove the contiguity, consider a simplex a in lZ a (P). We need to show that the vertices of a and h(a) 
span a simplex in 1Z a+2S (P). Clearly, all vertices of a are within a distance of each other. By definition 
of h, all vertices of h(a) are within distance a + 26. Let u be a vertex of a and h(v) be a vertex of /i(a) 
where v is a vertex of a. Then the Euclidean distance \\u — h(v) \\ is at most ||u — v\\ + \\v — /i(w) || < a + 5. 
Therefore, all vertices of a and h(o~) are within a + 25 distance. Hence, the simplex a and h(a) are faces 
of a common simplex in 1Z a+2S (Q) proving the claim of contiguity. ■ 

One may wonder how to efficiently construct the graph induced complexes in practice. Our experiments 
show that the following procedure runs quite efficiently in practice. It takes advantage of computing nearest 
neighbors within a range and, more importantly, computing cliques only in a sparsified graph. 

Let the ball B(q,5) in metric d be called the 5-cover for the point 9. A graph induced complex 
Q a (P, Q, d) where Q is a 5-sparse 5-sample can be built easily by identifying <5-covers with a rather standard 
iterative algorithm similar to the greedy (farthest point) iterative algorithm of |[T9l . Let Qi = {91, . . . , qi\ be 
the point set sampled so far from P. We maintain the invariants (i) Qi is <5-sparse and (ii) every point p € P 
that are in the union of 5-covers {J qe n- B(q, 5) have their closest point u(p) = argmin g q. d(p, 9) in Qi 
identified. To augment Qi to Qi+i = Qi U {qi + {\, we choose a point 9^+1 € P that is outside the 5-covers 
UqeQ B(q, 5). Certainly, 9j + i is at least 5 units away from all points in Qi thus satisfying the first invariant. 
For the second invariant, we check eveiy point p in the 5-cover of 9^+1 and update u(p) to be 9^+1 if its dis- 
tance to 9j + i is smaller than the distance d(p, v{p)). At the end, we obtain a sample Q C P whose 5-covers 
cover the entire point set P and thus is a 5-sample of (P, d) which is also 5-sparse. Next, we construct the 
simplices of Q a (P, Q, d). This needs identifying cliques in G a (P) that have vertices with different closest 
points in Q. We delete every edge pp' from G a (P) where v(p) = v(p'). Then, we determine every clique 
{pi, . . -Pk} in the remaining sparsified graph and include the simplex {v(pi), . . . , v{pk)} in G a (P, Q,d). 
The main saving here is that many cliques of the original graph are removed before it is processed for clique 
computation. We use the recently proposed simplex tree which computes cliques efficiently both time and 
space-wise (21. 

3 Surface point data 

In this section, we infer the geometry and topology of a surface through the graph induced complex. Let M 
be a smooth, compact surface embedded in IR 3 that has no boundary. We assume that M has positive reach 



p = p{M) which is the minimum distance of M to its medial axis. Let P be an e-sample of the metric 
space (M, dg) where dg is the Euclidean distance. Consider the graph induced complex Q a (P, Q, cLe)- In 
this section, the subset Q C Pis assumed to be a 5-sparse 5-sample of (P, de). 

Our result in this section is that under certain conditions on a, e and 5, Q a (P,Q,dE) captures the 
homology of M and contains the restricted Delaunay triangulation DelJM Q as denned below. The sparsity 
of Q turns out to be a crucial condition in the argument. 

Definition 3.1 Let Del Q denote the Delaunay triangulation of a point set Q C M 3 . The restricted Delaunay 
triangulation of Q with respect to a manifold M C M 3 , denoted Del|M Q> is defined to be the subcomplex 
of 'Del Q formed by all Delaunay simplices whose Voronoi duals intersect M. 

3.1 Topological inference from Q a (P, Q, d E ) 

Consider the sequence TZ a (P) -^ G a (P,Q,d) A K a+2S (P) in Proposition Ell When P, an e-sample 
of (M, ds), is sufficiently dense, it is known that i* : Hi(1Z a (P)) — > H\{TZ^{P)) is an isomorphism for 
appropriate a and j3. The homomorphism h* becomes injective if i* is an isomorphism since i* = j* o h*. 
If we can show that h* is also surjective, then h* becomes an isomorphism. We now show that h* is 
indeed surjective for Hi -homology and hence information about H\{M) can be obtained by computing 
V\\{G a {P, Q, dE))- First, we observe the following. Let P q C P be the points who have q G Q as the 
closest point. Notice that by the defmitino of h, h(P q ) = {q}. To prove that h* is surjective, it is sufficient 
to prove that the preimage of each 1-cycle in Q a {P, Q, d E ) contains a 1-cycle of lZ a (P). This, in turn, is 
true if the 1-skeleton of 7l a (P q ) is connected. 

Proposition 3.2 If the 1-skeleton oflZ a {P q ) is connected for all q £ Q, then h* is surjective. 

Proof: We show that the chain map h# induced by the simplicial map h is surjective. It follows that the 
homomorphism h* induced at the homology level is also surjective. Let c = qoqi+qiq2+- • -+QkQo be any 1- 
cycle in Q a (P, Q, dE)- The edges qi-iqi and qiqi+i have edges, say Pi-ip'i and Pip' i+ i respectively, in their 
preimage under h in TZ a {P). Consider a path 7$ between pi and p\ in 1Z a (P qi ) where h{pi) = h{p' i ) = qi. 
Such a path exists because lZ a (P q ) is connected for all q G Q. We have a 1-cycle 

c' = Pop'i + 7i + P1P2 + 72 + P2P3 -\ \- Ik + PkPo + 7o 

in 1Z a {P) so that h#{d) = c. This shows that h* is surjective in the first homology group. ■ 

The 1-skeleton of 1Z a (P q ) is connected if the union of balls B q = |J P B(p, t|) is connected because 
an edge P1P2 is in lZ a (P q ) if the respective balls B(pi,^) and B(p2,^) intersect. Let V q be the Voronoi 
cell of q in the Voronoi diagram Vor Q. Let M q = V q D M be the restricted Voronoi region. It turns out 
(we will prove it later in Proposition I3.41 i that if M q is contained in B q and M q is connected, then B q is 
connected. It may seem a priori that B q would contain M q if P is a dense sample. Unfortunately, that is not 
true as Figure |3]illustrates. To avoid such a case, we require that the Voronoi cells do not subtend very small 
angles between their facets which is ensured by the 5-sparsity of Q. Proposition [33] below uses 5-sparsity 
in a subtle way to prepare for the proof that B q contains M q . This result will also be used later to show that 
the graph induced complex Q a (P, Q, dE) in fact contains the restricted Delaunay triangulation Del|M Q- 

For a simplex a G Del|M Q, we call a ball B(c, r) a surface Delaunay ball of a if c € M and its 
boundary contains the vertices of a. 

Proposition 3.3 Let P be an e-sample of(M, dE), and Q a 5-sparse 5-sample of(P, dE)- Let a G Del|M Q 
be a restricted Delaunay triangle or edge with a vertex q G Q. Let c be the center of a surface Delaunay 
ball ofa.If8e<5< ^p(M), then there is a point p G P so that p G B(c, 4e) and q is the closest point to 
p among all points in Q. 




Figure 3: In a long thin Voronoi cell, B q may be disconnected and may not contain M q . 



Proof: See Appendix lAl for the proof. ■ 

Now, we are ready to prove that M q is contained in the union of balls |J p B(p,^). 

Proposition 3.4 If a > 12eand8e < 5 < -^p(M), then M q C U{ P eP } -^(p> § ) which implies that 
1Z a (P q ) is path connected if ' M q is path connected. 

Proof: Since P is an e-sample of M, Vx G M, there exists a point p G B(x,e) where p G P. Let 
P' q = P D (IJaeM B ( x i £ ))- Then, we have M q C U p6 p' ^(p> e ) for if x G M q> there exists p G P' q with 

p € B(x,e) requiring x G B(p,e). On the other hand, recall that P q = M q n P. Hence if p G P g \ P q , 
then S(p, e) contains some boundary point x G dM q . The point x belongs to a Voronoi facet in the Voronoi 
diagram of Q and hence B(x, \\q — x\\) is a surface Delaunay ball. By Proposition 13.31 we can find a 
point u G P q such that \\u — x\\ < 4e. Thus, B(p,e) C B(u, (4 + 2)e). Taking a > 12e, we get that 
M, C U PGP ^ B(p, e) C U p£Pq B(p, §). 

Since every ball in {B p \p G P g } intersects M q , we have that B q = U{ p eP } ^(p> § ) i s P at h connected 
if Mq is path connected. On the other hand lZ a {P q ) is path connected if B q is path connected proving the 
claim. ■ 

We can now present the main result of this subsection. 

Theorem 3.5 Let P be an e-sample of a smooth compact surface M embedded in H , and Q C P a 
5-sparse 5-sample of (P, fig?). For Vie < a < ^p and 8e < 5 < ^p, the map h* : h\\(TZ a (P)) — > 
\-\\{Q a {P, Q, dE)) is an isomorphism where h : lZ a {P) — > Q a (P, Q, cIe) is the simplicial map induced by 
the nearest point map v^ E : P — >■ Q. 

Proof: Since 5 < 0.18/0, we can assume each restricted Voronoi cell M q to be path connected lfl2l . This 
together with the lower bound on a imply that lZ a (P q ) is connected for each q G Q thanks to Proposition [ 
Consequently, Proposition I3.2l establishes that /i* is surjective. 



From Proposition 4.1 of |[T4l and its proof, we obtain the following: for any Ae < r < 2r < r' < 

Hi(ft r (P)) =* Hi(ft r '(P)) =* Hi(M) 



(1) 



where the first isomorphism is induced by the canonical inclusion i : 1Z T (P) ^-> 1Z r (P). Our assumption 



on the ranges of a and 5 implies the required conditions that 4e < a < |y§p and 4e < 5 < |a/| P- We 



claim that i* : Hi(7e a (P)) 
isomorphism. 



Hi(K a+2S (P)) induced by the inclusion i : TZ a {P) ->• K a+25 (P) is an 



First, note that this claim follows easily from Eqn £T|) if a < 5 by setting r = a and r' = a + 25(> 2r). 
Now assume that 5 < a. Consider the following sequence: 

K 5 {P) h K a {P) 4 K a+28 {P) ^ K ia (P). 

By Eqn £0), we have that the composition of inclusions ioi\ : 1Z S (P) — > 7Z a+2S (P) induces an isomorphism 
at the homology level. Hence i* is necessarily surjective. On the other hand, the composition of inclusions 
%2 o i : lZ a (P) — > 1Z 3a (P) induces an isomorphism at the homology level. Hence z* is necessarily injec- 
tive. Putting these two together, we have that i* is indeed an isomorphism. Therefore /i* is injective by 
Proposition I2.8I It then follows that /i* is an isomorphism as claimed. ■ 

Notice that the lower bound on 5 in Theorem I3.5l is not restricted by a. This means that one can have a 
dense input graph for a large a whose connectivity does not restrict the size of the subsample. 

In the next subsection, we show two examples of surface data where the graph induced complex has the 
correct Hi -homology with a considerably fewer simplices than the z^- witness complex, a modified witness 
complex suggested in (4l for capturing the topology correctly. 

3.2 Reconstruction of M using Q a (P,Q, d E ) 

In this subsection, we observe that the graph induced complexes can also be used for surface reconstruction. 
It is known that if P is dense and T is a simplicial complex with vertex set P which satisfies the following 
conditions, a simplicial manifold can be extracted from T that is homeomorphic to M EdU. The conditions 
are: (i) T is embedded in R 3 , (ii) all triangles in T have small circumradius with respect to reach and (iii) T 
contains the restricted Delaunay triangulation. We show that Q a (P, Q, d E ) contains the restricted Delaunay 
triangulation. We then prune Q a (P, Q,cIe) so that conditions (i)-(ii) are satisfied, but none of the restricted 
Delaunay triangles are deleted in the process which then ensures condition (iii). 

Theorem 3.6 For 8e < 5 < ^p and a > 8e, we have that Delji/ Q C Q a (P, Q,d E ) where P is an 
e-sample of (M, d E ) and Q C P is a 5-sparse 5-sample of (P, d E ). 

Proof: We will show that if 8e < 5 < ^p and a > 8e, then any triangle a € Del|M Q is in Q a (P,Q,d E ) ■ 
The theorem follows from this. 

Let a = <7i<72<73, and c the center of a surface Delaunay ball of a. By Proposition [331 there exists a point 
Pi G P in B(c, 4e) so that g, is the closest point in Q to pi for i = 1, 2, 3. It turns out that the interior of 
bounded cones used in the proof of Proposition I3.3l for q\, q2 and q^ are disjoint. Hence each point pi found 
in B(c, 4e) corresponding to qi is distinct from the other two. Therefore, if a > 8e, the vertices pi,p2 and 
P3 form a clique in G a (P) and hence the triangle q\q2qz is in Q a (P, Q, d E ). ■ 

The complex Q a (P,Q,d E ) may have intersecting triangles. We prune Q a (P,Q,d E ) to eliminate all 
such pairwise intersections while leaving the restricted Delaunay triangles in the complex. This ensures that 
the resulting complex embeds in M 3 and still contains the restricted Delaunay triangulation. Our simple 
observation is that if two intersecting triangles t\ and £2 do not intersect in a common face, one can decide 
locally which of the two can possibly be in a Delaunay triangulation. 

Observation 3.7 If V is the vertex set of two intersecting triangles t\ and £2 whose intersection is not a 
common face of both, then at least one of t\ and £2 is not in Del V. The triangle which is not in Del V 
cannot be in Del P where 7CP. 

One can check locally the Delaunay condition for t\ and £2 and decide to throw away at least one triangle 
which is not in Del V. This takes only constant time since V contains at most 6 vertices. Notice that no 
restricted Delaunay triangle can be thrown away by this process. After repeatedly pruning away one of 



the pairwise intersecting triangles, we arrive at a complex that embeds in ]R 3 and contains the restricted 
Delaunay triangulation Del|jvf Q. Next, we prune all triangles that have circumradius more than 26. Again, 
since the surface Delaunay ball of each restricted Delaunay triangle has circumradius at most 6 + e < 26, 
one is ensured that no restricted Delaunay triangle is eliminated. Assuming 6 to be sufficiently small, a 
sharp edge pruning and a walk on the outside of the resulting complex as described in ||3] [EH provides the 
reconstructed surface. The output surface enjoys one nice property that the triangles have bounded aspect 
ratios since they have circumradii of at most 26 and edge lengths of at least 6 (Q is 5-sparse). 

Theorem 3.8 If8s<6< ^fp, a > 8e, P is an e-sample of(M, cLe), and Q C P is a 6-sparse 6-sample 
of (P, cIe), then a triangulation T C Q a (P, Q,cIe) of M can be computed where each triangle in T has a 
bounded aspect ratio. 

Experimentally, we observe that surfaces can be reconstructed from a very sparse subsample with this 
strategy. Figure © presents two examples for surface reconstruction. The original sample P has 1, 575, 055 
points for the Fertility model and 1, 049, 892 points for Botijo model. The input graphs for the graph induced 
complex are constructed by connecting two points within distance of a = 0.45 for Fertility and a = 1.0 
for BOTIJO. The 2-skeleton of the Rips complex built on the input graph has 45, 788, 607 simplices for 
Fertility and 91, 264, 091 simplices for Botijo. The subsample Q consists of 3007 points for Fertility 
with 6 = 3.68, and 4659 points for BOTIJO with 6 = 4.0. The graph induced complex Q a (P,Q,dE) 
built on the subsample has : 3007 vertices, 9178 edges, 6304 triangles, 139 tetrahedra and no other higher 
dimension simplices for Fertility; 4659 vertices, 14709 edges, 10755 triangles, 718 tetrahedra, 5 4- 
dimensional simplices, and no other higher dimension simplices for BOTUO. The reconstructed surfaces 
from Q a (P,Q,dE) are shown in Figure @] For Fertility, it has 3007 vertices, 9039 edges and 6026 
triangles; for BOTIJO, it has 4659 vertices, 14001 edges and 9334 triangles. Evidently, the graph induced 
complex has only a few more simplices compared to the reconstructed surface. 

For a comparison, we also constructed the z/- witness complex suggested in [4] which also contains the 
restricted Delaunay triangulation Del|j\/ Q with v = (1, 6, 6, 4). The ^-witness complex for Fertility has 
3007 vertices, 35687 edges, 119237 triangles and 19874 tetrahedra; the v- witness complex for BOTUO has 
4659 vertices, 54648 edges, 180936 triangles and 29654 tetrahedra. The graph induced complex has much 
smaller size, but still captures f3\ (/?i = 8 for the Fertility, and /3i = 10 for the Botijo). 





(a) Fertility model (b) Botijo model 

Figure 4: Reconstructed surfaces for Fertility and Botijo models. 



4 Point data for more general domains 

In this section we consider domains beyond surfaces in R 3 . 

4.1 Manifolds 

Let M be a /c-manifold embedded in R n and P a discrete sample of (M, rig). We observe that the overall 
setup in section [3TT1 for inferring Hi-homology from the graph induced complex generalizes easily to higher 
dimensions. The inclusion map lZ a (P) '—t lZ a+25 (P) still induces an isomorphism at the homology level 
if a and S are chosen appropriately. In that case, the map h* : V\\{1Z a {P)) — > V\\(Q a {P, Q,d)) remains 
injective by the same argument as before. The main trouble arises when we try to prove that it is also 
surjective. Observe that, to prove the surjectivity of h*, we used the fact that the restricted Voronoi cell 
M q = Vq n M in Vor Q is connected (Proposition I3.4I ). Unfortunately, this is not necessarily true in high 
dimensions given the counterexamples in [@][T0]]. To overcome this impediment we change the distance 
function replacing the Euclidean distance with the graph distance while building Q a (P, Q,d). Specifically, 
we still consider G a (P) to be the graph connecting points in P with Euclidean distance a or less, but take 
Q to be a 5-sparse <5-sample of (P, da) where the graph distance dg = dQa^ is defined with the Euclidean 
lengths as the edge weights. Then, we consider Q a (P, Q, do)- 

As before, let P q C P be the set of points nearest to a point q € Q with respect to da- The modification 
in distance function immediately allows us to claim that the 1-skeleton of lZ a (P q ) is connected, which was 
needed to claim that h* is surjective. 

Proposition 4.1 !Z a {P q ) is connected, and thus h* is surjective. 

Theorem 4.2 Let P be an e-sample of an embedded smooth and compact manifold M with reach p, and 
Q a 5-sample of(P,d G ). For Ae < a, 6 < \J\fh the map h* : H x {n a {P)) -> H 1 (g a (P,Q,d G )) is an 
isomorphism where h : lZ a (P) — > Q a (P,Q,dQ) is the simplicial map induced by the nearest point map 
v dG : P -> Q. 

4.2 A leaner subsampling for Hi 

In this subsection we show that the subsample Q can be made leaner. The main insight is that we can 
define a feature size larger than the reach which permits us to subsample more sparsely with respect to 
this larger feature size. Gao et al. liTTl considered a similar feature size for the same reason of requiring 
sparser sampling for a two dimensional shape. Here we show that such a sparser sample is also adequate for 
determining Hi of manifolds in high dimensions. Our experimental results in Figure |2] suggest that one can 
obtain information about Hi from a very sparse sample in practice. 

Let K, be a simplicial complex with non-negative weights on its edges. We define homological loop 
feature size as 

1 inf{|c|, c is non null-homologous 1-cyclein/C} 



hlfe(/C) = { 2 .. 

[ oo it no such c exists. 

This feature size is very similar to the systolic feature size sfs(X, d) of a compact metric space (X, d) ifTTl 
which is the length of the shortest non-contractible loop in X. Our definition of hlfs when applied to a 
metric space (X, d) becomes larger than or equal to sfs(X, d) . Notice that every loop of K with length less 
than 2hlfs(/C) is null-homologous in K. 

Let Q C P be a 5-sample of (P, d G ) as before but with 5 < ihlfs(^ Q (P)) - \a. Let h : K a (P) -> 
G a {P,Q,dc) be the simplicial map as defined earlier. We aim to show that the induced homomoiphism 
h* on the first homology is injective. Since we use graph distances, Proposition 14.11 remains valid and 
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hence h* remains to be surjective. However, we cannot claim z* : Hi(JZ a (P)) — > Hi(JZ a+2S (P)) to be an 
isomorphism because 5 could be larger than required. Thus, we cannot use i* to infer h* to be injective as 
before. Nevertheless, we can prove the following result using a different approach. 



Theorem 4.3 // Q is a 5-sample of (P, d G ) for 5 < ±Ms(K a (P)) 
Hi(^ a (-P, Q,dc)) is an isomorphism. 



\a, then K : Hi(K a (P)) 



luv 





(a) 



7fc 

v = Vk 
(b) 




Figure 5: (a) j uv makes a cycle with tt(u, a), ir(v, b) and ab, (b) j uv as a sum of unicolored chains and 
bicolored edges, (c) converting j uv (shown dotted) to j uv , (d) a diamond of j uv . 



Proof: We only need to show that h* is injective, as its surjectivity follows from Proposition 14.11 To show 
the injectivity, it suffices to show that h* has a trivial kernel. Let a be any triangle in Q a (P, Q, do)- If under 
the chain map h#, every cycle in the preimage of the boundary cycle da is null homologous, then every null 
homologous cycle in Q a {P, Q,dc) has only null homologous cycles in its preimage. This is true due to the 
fact that a bounded cycle is a sum of boundaries of triangles, and that the chain map h# is surjective (see 
the proof of Proposition I3.2I ). Below we show that under the chain map h#, every cycle in the preimage of 
the boundary cycle of any triangle indeed is null homologous. It would then follow that the kernel of h* is 
trivial. 

Let 7 be any cycle in the preimage of da. We have 7 € Y^uveda h~Z (uv) where uv be any edge of of a. 
Let ~f uv be any maximal subpath of 7 so that h#(^ uv ) = uv (Figure 13 a)). For each such j uv , we construct a 
cycle £ uv so that £ uv is null homologous and 7 is homologous to ^ £ uv . Therefore, our problem of showing 
7 is null homologous reduces to the problem of showing every £ uv is null-homologous. 

We construct £ uv as follows. By the construction of Q a {P,Q,dc), there is a triangle abc G lZ a (P) 
such that h(abc) = a with h{a) = u, h(b) = v. Consider the shortest paths n(u, a) and ir(v, b) in G a (P) 
from u to a and from v to b respectively. Observe that all vertices in ir{u, a) and 7r(u, b) are mapped to 
u and v respectively by h since we are using the graph-induced distance da to construct Q a (P,Q,dG)- 
Take £ uv to be the chain n(u, a) + ab + tt(v, b) + j uv ; refer to Figure |5ja). With this choice, we have 
7 = Y^ (-uv + d(abc) and hence 7 is homologous to ^ £ uv as promised. To prove £ uv null-homologous, we 
construct a homologous path 7 U1) to j uv which gives a homologous cycle £ uv to £ uv . We then prove that £ uv 
is null-homologous. 

Call an edge e = (x,y) in TZ a (P) bicolored if its two end-points are mapped to two distinct vertices 
by h; otherwise, e is unicolored. A 1-chain from 1Z a (P) is unicolored if it has only unicolored edges. For 
simplicity, we assume that vertices from r y uv are all contained in h~ l (u) U h~ l {v) because ^ uv is always 
homologous to a path containing vertices mapped only to u or v. In this case, j uv can be decomposed into 
a set of bicolored edges {yoxi, 2/1X2, . . . , yk-iXk} together with a set of unicolored chains {70, 71, • • • , 7fc} 
such that dji = Xi + yi. In particular, ^ uv can be written as 



7« 



7o + 2/0^1 + 71 + yix 2 H h 2/fc-iXfc + Ik- 



(2) 
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See Figure [5jb) for an illustration. 

Claim 4.4 Let ji be a unicolored chain with two boundary points Xi , yi so that for any simplex r € ji, 
h(r) = u. Then, 7, is homologous to the chain 7$ = ir(xi,u) + ir(u, y,). 

Given the chain ^y uv with subchains 7$ as in equation |2 we convert it to a homologous chain 

%v = 70 + Voxi + 71 H h g/fc-iXfc + Ik- 

Replace 7 m , with •y uv in £ TO to obtain a homologous cycle £ uv ; refer to Figure [5] Observe that £ uv is the 
sum of cycles (diamonds) that have two unicolored chains and two bicolored edges as shown in Figure Eld). 
Such a cycle c has length at most 45 + 2a. This is because each unicolored chain in c has at most two 
shortest paths of the form ir(u, X{) and ir(u, yi) (or ir(v, xi) and n(v, yi)) that have lengths 25 or less (Q is 
a 5-sample of (P, da)), and the two bicolored edges have lengths at most 2a in G a (P). The cycle c is null 
homologous because its length is 

\c\ <45 + 2a< 2h\k(TZ a (P)), given that 5 < ihlfs(ft Q (P)) - \-a. 

It follows that £ uv is null homologous as we are required to prove. 

We only need to show Claim |4~41 to finish the proof. Let cc, = Po,Pi, • • • ,Pm = y-i be the sequence 
of vertices on the unicolored chain (path) 7$. Consider the shortest paths ir(pi,u) for each pi on this path. 
The length of the cycle Zi = ir(u,pi) + PiPi+i + n(u,Pi+i) is at most 25 + a for each i € [0, m — 1]. 
Therefore, it is null homologous by our assumption. We have 7$ + 7$ = ^I^-o Zi = *-*■ Therefore, 7$ and 7j 
are homologous. ■ 

4.3 Point data for compact sets 

So far we have focused on Hi -homology. In this section we extend the domain to compact subspaces of 
Euclidean spaces and consider homology groups of all dimensions. This generality comes at the expense of 
additional computations. Unlike previous approaches that allow us to infer the Hi -homology of the sampled 
manifold by computing directly the same for the graph induced complexes, now we need to compute the 
persistent homology |[T5l induced by simplicial maps. The well-known algorithms for computing persistent 
homology lfI31 work for maps induced by inclusions. In a contemporary paper |[T3l . we present an algorithm 
that can compute the persistent homology induced by simplicial maps. 

Let 1 C 8" be a compact set and X x be its offset with A > 0. Since it is difficult to compute Hfe(X) 
from a sample J9|, we aim for computing the homology groups Hfc(X A ) for the offset X x . Let wfs(X) 
denote the the weak feature size which is defined as the smallest positive critical value of the distance 
function to X (H. We prove that the persistent homology of the graph induced complex defined with either 
Euclidean or graph distance d provides the correct homology of X x where < A < wfs(X). Specifically, 
the image of h* : H k (g a (P, Q, d)) -> H k (Q a ' ' (P, Q' , d)) induced by a simplicial map h : G a (P, Q, d) -> 
Q a (P, Q', d) becomes isomorphic to Hfc(X A ) for appropriate a and a'. We recall the following result from 
0- 

Proposition 4.5 If the sequence of homomorphisms A^-B^-C^-D^-E^-F between finite 
dimensional vector spaces satisfies that r&nk(A — > F) = rank(C — > D), then rank(B — > E) = rank(C — > 
D). 
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Let Q and Q' be subsamples of P where Q is a <5-sparse <5-sample and Q' is a ^'-sparse 5 -sample for 
5 > 5. Consider the interleaving sequence between the graph induced and Rips complexes, 



TZ a (P) ( - > Tl a+25 (P) c n > 7£ 4 ( Q + 2<5 ) (P) c n > 7e 4 (a+ 25 )+25 (p) q) 




£ a (P, Q, d) h -> S 4 ("+ 25 )(P, Q', d) 



where i\, 12,13, ji and j'2 are inclusions and h = h^ov^o j\. By Proposition 12.81 h\ and /12 are simplicial 
maps . Therefore, h is also a simplicial map as composition of simplicial maps. In particular, h is the 
simplicial map induced by the vertex map that maps each point q € Q to its closest point q' G Q' in Q'. We 
prove that im h* = H fc (X A ) where h* : H fc (£ a (P, Q, d)) ->■ H fc (0 4 ( a+M ) (P, Q', d)) and e , a and 5 fall in 
appropriate ranges. 

Theorem 4.6 Lef X C R n Z>e a compact space. Let < e < g wfs(X) and P be an e-sample of (X, d#). 
Let Q be a 5-sparse 5-sample of (P, d) and Q' be a 8' -sparse 5 -sample of (P, d) where d is either Euclidean 
or graph distance and 5 > 5. 

If 2s <a< j(wh(X) - e) and (a + 25) + \# < |(wfs(X) - e), then im h* = H k (X x ) (0 < A < 
wfs(X)) where h* : H k {G a {P, Q, d)) -> H k (Q 4( - a+2 ^ (P, Q', d)) is induced by h in diagram\3\ 

Proof: The diagram [3] is not commutative in general. However, it is commutative at the homology level. 
Proposition I2.8l makes the two triangles at the left and right commutative. The middle square commutes by 
definition of h. Now consider the sequence, 

H k (K a (P)) -^ H k (g a (P, Q, d)) ^-^ H k (K a+25 (P)) (4) 

^— H k {TZ^ a+26 \P)) — H fc (g 4 ( Q+2<5 )(P,Q',d)) -^ H k (TZ^ a+25 ^ +28 ' (P)) 

Consider the sequence of inclusions at the upper level of the diagram (1). Since a > 2eand (a + 25) + ^5 < 
|(wfs(X) — e), we have that 

im (j 3 o i 2 o »!), =* H fc (X A ) and im (i 3 ), =* H fc (X A ) 

by Theorem 3.6 of O. Considering the diagram in (1) and the sequence in dU we have 

H fe (X A ) = im (i 3# o z 2# o i u ) ^ im((j 2 * ° ^2*) ° «2* ° (ji* ° hi*)) - im «2* (5) 

Letting /i = h 2 i 2 o j 1; the rightmost isomorphism in © allows us to claim that im h* = H k (X x ) by 
applying Proposition !4.5l to the sequence ©. ■ 

5 Conclusions 

In this work, we investigated the graph induced complex that can be built upon a given point cloud data 
and a suitable graph connecting them. This complex, to some extent, has the advantages of both Rips and 
witness complexes. We have identified several of its topological properties that can evidently be useful in 
extracting information from point data even in high dimensions. 
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In section I4.3L we have shown how one can infer the homology groups of a compact set using the 
persistent homology of a pair of graph induced complexes constructed with two values of 5. One can 
consider a filtration of Q a (P,Q,d) with Q being sparsified for increasing values of a and 6. Then, one 
can obtain a persistence diagram ifTTI out of this 'full filtration' using our recently proposed algorithm for 
computing the topological persistence for filtrations connected with simplicial maps |[T3ll . The algorithm will 
collapse vertices progressing through the nitration and hence will keep the size of the complex in question 
contained. Relating this persistence diagram to that of a filtration obtained by a related Rips filtration is an 
interesting question. In a subsequent work, we have addressed this question in lfT3l . 

Finding other applications where the graph induced complex becomes useful also remains open for 
further investigations. 
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A Proof of Proposition 133 



First, we present an elementary geometric result that we need to use in the proof. Let C(o, v, a) denote the 
cone with apex o, axis in the direction of v and aperture 2a. 

Claim A.l Given a ball B with radius r and center c, let q be an arbitrary point on the surface of sphere. 
Consider the two nested cones C\ = C(c, cq, a) and C2 = C(c, cq, 2a) with the same axis. Let p be any 
point from the intersection of the ball and the inner cone; that is, p 6 B n Ci. Let x be an arbitrary 
point from the boundary of B outside the outer cone; that is, x £ C2, and x € dB. Then we have that 

\\p — q\\ < \\p ~ x \\- 

Proof: Denote a q := Zpcq and a x := Zpcx. Because x is outside of the outer-cone with aperture 4a, and 
p is inside of the inner cone of aperture 2a, we have that we have a q < a < a x . Now consider the triangle 
Apcq. By the Cosine Law, we have that 

\\p ~ q\\ = \\p ~ c \\ +ll c — q\\ ~ 2||p — c\\ ■ \\c — q\\ cos(ad) = \\p — c\\ + r — 2r\\p — c\\ cos(a g ). 
Similarly, consider the Apox, and we have 

\\p — x\\ = \\p — c\\ + \\c — x\\ — 2\\p — c\\ ■ \\c — x\\ cos(a p ) = \\p — c\\ +r — 2r\\p — c\\cos(a x ). 

Since < a q < a x < n, we have \\p — q\\ < \\p — x\\. ■ 

Now consider the surface Delaunay ball B c = B(c, r) that passing through the vertices of the simplex 
a and containing no other points from Q. Recall that q is an arbitrary vertex of a. Since all other vertices 
of a are at least 5-Euclidean distance away from q, we then have that the intersection of B c with the cone 
C(c, cq, 2 arcsin j-) contains no point from Q other than q. By applying Claim lATTI with a = arcsin ^-, we 
then obtain that: 

Corollary A.2 If there exists a point p € P such that p € B c n C(c, cq, arcsin ^) n M, then q must be the 
closest point to p among all points in Q. 

In what follows, we will show that a point p € P satisfying the conditions in Corollary |A.2| as well as 
that p e B(c, 4e) indeed exists when 8e < 5 < ^p(M). This will then prove the proposition. Specifically, 
we will first identify a sample point p£P, and then we will show that p satisfies the requirements of the 
proposition. 

Identifying a point p € P. Let B = B(m,p) and B / = B(m',p) be two 
balls tangent to M at c; assume without loss of generality that B is inside of 
M and B < is outside. Locally around c, the surface M is sandwiched between 
B and B i. Now consider the plane V = span{o, o',q}; note that c also lies 
in V. Denote B c<v = B c n V, B oP = B n V and B , jV = B , n V. Let x 
be the intersection point of B c p and B p that is on the same side of the line 
passing through 00' as the point q. Similarly, let y be the intersection point of 
B c -p with B i-p on the same side of the line 00' as q. Obviously, q lines on the 
arc xy that avoids B -p and B < -p. See the right figure for an illustration where 
the shaded region is B C: p. Set 9 := Zxcy; easy to see that 9 = Zxoc = Zyo'c. Hence we have that 

9 ||c— a:|| r 

S111 2 — 2p ~ Tp- 

Now consider the segment cq and the point w € cq such that the length of cw is a value £ which we 
will set later. As w is contained in the cone with apex c and aperture 9 in the plane V, the ball B(w,£ sin 9) 
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will intersect both segment ex and cy, thus intersecting both B and B >. Since B is inside the surface M 
and B , outside, it follows that B(w,£sin9) nM/l Pick any point p' G B(w,£s'm9) n M. By the 
e-sampling condition of P, there must exist a sample point p € P such that ||p — p' \\ < e. In other words, 
there is a sample point p G P such that p G -B(io, £ sin + e). 

The requirements on p. We now need to show that the parameter £ can be chosen such that the point p 
satisfies all the requirements from the proposition. In particular, we need the following: 

C-l pG B(c,4e); and 

C-2 q is the closest point to p among all points in Q. 

Now set r = arcsin - sin f +e . Obviously, the ball B(w,£ sin 9 + e) (and thus the point p) is con- 
tained inside the cone C(c, cq, r). Observe that by Corollary IA.2I condition (C-2) is satisfied if (C-2.a) 
B(w,£ sin 9 + e) G B(c, r) (implying that p G M), and (C-2.b) r < arcsin ^ (implying that p is contained 
in the inner cone C(c,cq, arcsin ^-)). 

The existence of a valid L What remains is to find a value for £ so that (C-l), (C-2.a), and (C-2.b) are 
all satisfied simultaneously. Note that since ||io — c|| = £, we have that \\p — c\\ < £ + £sin9 + e. Hence 
condition (C-l) is satisfied if £ + £sin9 + e < 4e. Since sin^ < 2sin| = -, (C-l) holds as long as the 
following inequality holds. 

Since 5 > 8e, if (C-l) holds, then we have that ||p — c\\ < 4e < | < r, which implies (C-2. a). Now 

consider condition (C-2.b), which holds if ^ sin / +£ < |:. Since 5/2 < r < 5 + e and 8e < S < 2p/27, we 

have that: 

r 5 + e ^ 5 + 6/8 1 6 ^ 6 
- < < — < - < — < — • 

p p p 4 AS 2r 

That is, ^ - sin 9 > ^ - - > (recall that sin 9 < 2 sin 9 = r/p). Hence condition (C-2.b) holds if 

l>^r(> , £ . )■ (7) 

77 TT 8111 ft 

Putting Eqns © and (O together, we have that as long as the value £ satisfying the following inequality: 

J> r — — 1_|_1 

2r p ^ P 

then conditions (C-l) and (C-2) will be satisfied, and there exists a point pgPas stated in the proposition. 
Given that 8e < S < 2p/27, we can show that valid £ exists. For example, for £ = ^, inequality in Eqn 
© holds as 



e e e e 36e 

< —l ^ < oT< 



and 



6 r - 8 5+e — AX lg A 2.2£ "' 1Q 



3e 3e 3e 3e 36e 

> " TTT > OT > 
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